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Abstract. Compactness of some stars is explained if they 
are strange stars (SS) as shown by Dey et al. (1998) (D98) 
and Li et al. (1999a). One of these compact star candi- 
dates is the SAX J1808.4-3658 (SAX in short) believed 
to be an important link in the genesis of radio pulsars. 
SS have also been suggested for bursting X-ray pulsars 
(GRO J1744-28, Cheng et al. 1998), from quasi-periodic 
oscillations (QPO) of X-ray binaries (4U 1728-34, Li et al. 
1999b) and from peculiarity of properties of radio pulsars 
(PSR 0943+10, Xu et al.1999; Kapoor et al. 2000). We 
now extend the calculation to include high temperatures 
upto T = 70 MeV - 8 x 10 11 °K and hnd that the nature 
of the mass (M) and radius (R), derived from astrophys- 
ical data, is still retained. The entropy is calculated and 
matches onto that calculated from hadronic models thus 
supporting the idea that the quark-hadron transition may 
be continuous. 



to these stars which emit X-rays, generated presumably 
due to accretion from their binary partner. 

During the genesis of these stars higher T may be en- 
countered and in this paper we deal with a generalized 
case of an object at a uniform T. We show at upto T = 
70 MeV a self sustained system can be supported by the 
parameter set eosl (D98) (called SSI in Li et al. 1999a). 

The changes in star masses are shown in the range of 
T mentioned above. The conclusions of D98 and Li et al. 
(1999a, 1999b) are still valid for the finite T cases. Sec- 
ondly, parameters of the single particle potential at finite 
T are tabulated. In particular the entropy is studied and 
this can be compared with hadronic models. The com- 
parison shows that the entropy may indeed be continuous 
supporting the idea of a continuous phase transition be- 
tween hadrons and quarks. 



2. Strange stars at finite temperatures 
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1. Introduction 

A calculation for cold strange stars (D98) enabled us to 
draw conclusions about chiral symmetry restoration in 
QCD when the EOS was used to get SS fitting definite 
mass-radius (MR) relations (D98; Li et al. 1999a; Li et 
al. 1999b; Ray et al. 2000b) . The empirical MR relations 
were derived from astrophysical observations like luminos- 
ity variation and some properties of quasi-periodic oscilla- 
tions from compact stars. The calculations are compared 



The interesting point made in D98 is that starting from 
an empirical form for the density dependent masses of the 
up (u), down (d) and strange (s) q-s given below one can 
constrain the parametric form of this mass from recent 
astronomical data[]: 

Mj = mi + Mq sech(^ — ), i = u, d, s. (1) 
Po 

Restoration of chiral quark masses at high density is incor- 
porated in this model. Using the model parameter (v) for 
this restoration one can calculate the density dependence 

1 pB = (pu + Pd + Ps)/3 is the baryon number density, po — 
0.17 fm -3 is the normal nuclear matter density, and v is a 
numerical parameter. The current q masses mi used in the 
following are 4, 7 and 150 MeV for u, d and s respectively. 
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of the strong coupling constant (Ray et al. 2000b). In other 
words the masses of stars in units of solar mass, (M/Mq), 
found as a function of the star radius i?, calculated using 
the above Eq. |], - produces constraints which enable us to 
restrict the parameter v. At high pb the q- mass Mj falls 
from its constituent value Mq to its current one mi. The 
parameter Mq is taken to be 310 MeV to match up with 
constituent q- masses assuming the known fact that the 
hadrons have very little potential energy. The results are 
not very sensitive in so far as changing Mq to 320 MeV 
changes the maximum mass of the star from 1.43735 Mq 
to 1.43738 Mq and the corresponding radius changes from 
7.0553 kms to 7.0558 kms. 
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Fig. 1. The smooth restoration of chiral symmetry inside 
the star for each of the u, d and s flavours (note that this 
is for the zero temperature result). 

It is interesting to plot the up (u), down (d) and strange 
(s) q- masses at various radii in a star. This is done with 
strong coupling constant ao — 0.2, chiral symmetry restora- 
tion parameter v = 1/3 and the QCD scale parameter 
A = 100 already discussed in D98 & Li et al. (1999a, 
1999b). Fig. |l| shows that the quarks do not have the con- 
stituent masses as in zero density hadrons nor do they 
have the current masses of the bag model. Upto a radius 
about 2 kms the quarks have their chiral current mass but 
in the major portion of the star their masses are substan- 
tially higher. At the surface the strange q- mass is about 
278 MeV and the u, d q- masses - 130 MeV. 

We use the large Ac (colour) approximation of 't Hooft 
(1974) for quarks, where quark loops are suppressed by 
1/A C and the calculation can be performed at the tree level 
with a mean field derived from a qq interaction (Witten 



Fig. 2. The equation of states for different T. It is to be 
noted that beyond T = 70 MeV, the zero of the pressure 
cannot be attained. 

1979). This was done for baryons (Dey et al. 1986; Dey et 
al. 1991; Ray et al. 2000a) and extended to dense systems 
like stars (D98). Following is the Hamiltonian, with a two- 
body potential Vy : 

H = 5>, Pl + (3 tMl ) + £ (2) 

% i<j 

The vector potential in Eq. ^ between quarks originate 
from gluon exchanges, and the A-s are the color SU(3) ma- 
trices for the two interacting quarks. In the absence of an 
accurate evaluation of the potential (e.g. from large A c 
planar diagrams) we borrow it from meson phenomenol- 
ogy, namely the Richardson potential (Richardson 1979). 
The potential reproduces heavy as well as light meson 
spectra (Crater et al. 1984). It has been well tested for 
baryons in Fock calculations (Dey et al. 1986; Dey et 
al. 1991). Recently, using the Vlasov approach with the 
Richardson potential, Bonasera (1999) finds a transition 
from nuclear to quark matter at densities 5 times pQ. 

We had to calculate the potential energy (PE) contri- 
bution in two steps : the single particle potential, Ui{k) 1 for 
momentum k, is first calculated and this is subsequently 
integrated to get the PE. The Ui(k) is needed for doing 
finite T calculation. 

The Ui(k) is parametrized as sum of exponentials in k 
(i.e., Ui(k) = — exp(ao + a\k + 02k 2 ) for a given flavour), 
where the parameters for a given set, reported in this pa- 
per, with v = 1/3, A = 100 MeV, a = 0.2 and a typical 
density 10po are given in Table |. 
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Table 1. Parameters of Ui(k) = —ex-p(ao i + ai i ki + a2ikf) 



i 


ao 


ai 


«2 


u 


-2.3524 


0.10162 


-0.09832 


d 


-2.3486 


0.10281 


-0.10003 


s 


-2.3109 


0.12042 


-0.11023 



Finite temperature T(= 1/(3) can be incorporated in 
the system through the Fermi function : 

FM{k, T) = =-— r= (3) 

exp[/3(e - e F )} + 1 

with the flavour dependent single particle energy 

a = ^k 2 + M i ( P ) 2 + U i {k). 

Now we evaluate 



(4) 



7 



(e)k 2 FM{k, T)dk 



For cj)(e) = 1 we get the number density and for <j){e) = 
e, the energy density. 7 = 6 is the spin-colour degeneracy. 

Even at very high T which is around 70 MeV, the 
chemical potential is very large, of the order several hun- 
dred MeV. The entropy is calculated as follows: 



s{T) 



— / k 2 \FM(k, T)ln(FM(£;, T)) 
tH Jo 



+ (1 - FM(k, T))ln(l - FM(k, T))]dk (5) 
The pressure (P) is calculated from the free energy 

f = e-Ts 

as follows: 



We find P = points only upto T = 70 MeV on plotting 
P as a function of p at different T (Fig. 0) . 



3. Results and discussions 

With this EOS at finite T, the TOV equation is solved to 
find out the MR relationship of the strange star at finite 
T. Here we find that the stellar mass and radius decrease 
with increasing T and the central energy density goes upto 
22 times the normal saturation density of nuclear matter. 
No star can be formed beyond T = 70 MeV since the star 
surface must have pressure zero for SS. The MR curves 
for the stars at the three values of T are given in Fig. || 
showing that the radii and the masses decrease slightly 



with T. The conclusions of D98 and Li et al. (1999a & 
1999b) remain unchanged so far as compactness of the 
stars is concerned. 

We thank the anonymous referee for pointing out to 
us that the decrease in the maximum mass and the corre- 
sponding radius i? max of the star at finite T is somewhat 
surprising. We offer explanations of this fact as follows : 

(1) We show in the Fig. ^ that the mass of the star 
with a given radius increases with T. This is clear from 
the vertical line which connects the maximum mass for T 
= 70 MeV with i? max = 6.076 km to the star mass for T 
= for the same radius. This is in accord with the usual 
expectation that the star mass should increase with T. 

(2) With increasing T, the size of any self sustaining 
system decreases due to the restriction placed on the en- 
ergy balance by the increased thermal energy. The shrink- 
ing of a self-sustained system with increasing T is also 
seen for the Skyrmion (Dey and Eisenberg 1994). Inter- 
estingly the latter is also an example of the success of the 
large N c phenomenon. Note that P in our system (Fig. ||), 
or the Skyrmion, is calculated self consistently whereas in 
the earlier literature on strange stars at finite T (Kettner 
et al. 1995) employing bag model, the variation of the bag 
pressure with T is neglected. It is pointed out in (Chmaj 
and Slominski 1989) at low P the bag pressure dominates 
and as P grows the results are close to the free relativis- 
tic gas limit. With an unchanged bag parameter therefore 
bag model calculations find an almost unchanged mass 
and radius for the stars at finite T . 

(3) As thermal energy increases the binding energy per 
baryon decreases. This supports the argument given in (1) 
and (2) above that as T decreases the masses correspond- 
ing to the same radius will decrease. 

The entropy S increases with T and decreases with 
density. Comparable to our entropy is that calculated by 
Das, Tripathi and Cugnon (1986) (DTC) for interacting 
hadrons. The results of DTC checks with experiment. The 
minimum T considered by DTC is 10 MeV and S/A is a 
little less than 1 at a density 2.5 times po 

In Fig. U we have plotted S/A as a function of the star 
radius for T = 20, 50 and 70 MeV. It is interesting to 
see here that entropy is maximum at the surface showing 
that the surface is more disordered than the core. The 
extrapolated entropies from the plots of DTC at 5 po agree 
with our S/A at the stellar surface for all the three cases. 

Recently Glendenning (Glendenning 2000) has argued 
that the SAX could be explained as a neutron star rather 
than bare SS, not with any of the existing known EOS, but 
with one based only on well-accepted principles and hav- 
ing a core density about 26 po- Of course, such high density 
cores imply hybrid strange stars, subject to Glendenning's 
assumption that such stars can exist with matching EOS 
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having minimum at about 4.8 po whereas in Glendenning 
(2000) one of the assumptions is that strange matter has 
no bound state. 




In summary, we find that beyond T = 70 MeV, the 
EOS has no zero pressure point. A self-bound star cannot 
exist in our model at higher T. The entropies at T = 20, 
50 and 70 MeV (intermediate values are quite obvious) 
match onto hadronic entropies at corresponding T, sug- 
gesting the possibility of smooth phase transition between 
the hadronic and the quark states. 
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Fig. 3. MR curves for T = 0, 20, 50 & 70 MeV. Core 
densities for T ^ are 16.5, 19.5 and 22 po- We have 
drawn a vertical line at i? max for T = 70 MeV to show 
that the mass of star increases as T increases, for fixed 
radius. 
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Fig. 4. The entropy per quark at different T(20, 50 & 70 
MeV) as a function of the star radius. 



for two phases. There is the further constraint that if the 
most compact hybrid star has a given mass, all lighter 
stars must be larger. It was found in Li et al. (1999b) 
that the star 4U 1728—34 may have a mass less than that 
of SAX and yet have a radius less than -Rsax- Another 
serious difference is that in our model the strange quark 
matter EOS derived, using the formalism of large N c ap- 
proximation, indeed shows a bound state in the sense of 
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